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Summary 


Molecular changes elicited by plants in response to fungal attack and how this affects plant- 
pathogen interaction, including susceptibility or resistance, remain elusive. We studied the 
dynamics in root metabolism during compatible and incompatible interactions between chickpea 
and Fusarium oxysporum f. sp. ciceri (Foc), using quantitative label-free proteomics and NMR- 
based metabolomics. Results demonstrated differential expression of proteins and metabolites 
upon Foc inoculations in the resistant plants compared with the susceptible ones. Additionally, 
expression analysis of candidate genes supported the proteomic and metabolic variations in the 
chickpea roots upon Foc inoculation. In particular, we found that the resistant plants revealed 


significant increase in the carbon and nitrogen metabolism; generation of reactive oxygen 
species (ROS), lignification and phytoalexins. The levels of some of the pathogenesis-related 
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Introduction 


Chickpea (Cicer arietinum L.) is the second most widely grown 
legume in the world. As a top producer, India contributes about 
90% of global chickpea production (http://faostat.fao.org/site/ 
339/default.aspx). Chickpea is mainly used as a primary vegetar- 
ian source of human dietary protein and, thus, is of significance to 
food and nutritional security in the developing world. However, 
due to widespread occurrence of fungal pathogens, such as 
Fusarium oxysporum and Ascochyta rabei, the yield of chickpea 
has been constrained in spite of successive efforts of national and 
international breeding programmes. Annual yield losses due to 
wilt disease alone have been estimated to range from 10 to 90% 
(Anjaiah et al., 2003; Jimenez-Diaz et al., 1989). Wilt is caused by 
eight races of Fusarium oxysporum f. sp. ciceri (Foc) affecting all 
the major chickpea growing areas (Gurjar et al., 2009). Foc 
infects roots and clogs the xylem, resulting in the obstruction of 
nutrient supply as wilting progresses, ultimately leading to plant 
death. This fungus can survive for many years in soil even without 


proteins were significantly higher upon Foc inoculation in the resistant plant. Interestingly, results 
also exhibited the crucial role of altered Yang cycle, which contributed in different methylation 
reactions and unfolded protein response in the chickpea roots against Foc. Overall, the observed 
modulations in the metabolic flux as outcome of several orchestrated molecular events are 
determinant of plant’s role in chickpea—Foc interactions. 


its host and, hence, poses a serious challenge for disease 
management (Haware et al., 1996). 

Many studies have been carried out to identify the molecular 
basis of Foc resistance or susceptibility in chickpea using various 
approaches including gene mapping, candidate gene identifica- 
tion, differential expression and biochemical analysis postfungal 
infection (Ashraf et al., 2009; Giri et al., 1998; Gowda et al., 
2009; Gupta et al., 2010, 2013a; Gurjar et al., 2012; Nimbalkar 
et al., 2006). However, genomic scale dynamics of plant-fungus 
interaction is poorly understood. Further, due to variations in 
chromosomal rearrangements between nonmodel and model 
plants, differences in the plant-pathogen interaction with respect 
to signalling and immunity events are also expected and evinced 
(Gupta et al., 2013a). Therefore, studies on nonmodel plants 
using unbiased modern high-throughput technologies are 
required to improve our knowledge of the plant-fungus interac- 
tions (Kushalappa and Gunnaiah, 2013; Mehta et al., 2008). In 
recent years, a combination of metabolomics and gene-expres- 
sion analysis has been employed to understand such interactions 
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(Liu et al., 2010) and the effects of gene manipulation on the 
systems metabolic changes of Fusarium graminearum (Chen 
et al., 2011). A few studies using functional genomics (Cho et al., 
2012; Golkari et al., 2007), proteomics (Lee et al., 2009; Yang 
et al., 2010) and metabolomics (Bollina et al., 2011; Kumaras- 
wamy et al., 2011) approaches have also been reported on plant- 
pathogen systems of Fusarium head-blight and cereal crops. 

As plant—fungus interactions are complex in nature and factors 
leading to resistance or susceptibility in plant remain largely 
obscure, the objective of current study was to assess overall 
modulations in the levels of proteins and metabolites in chickpea 
roots upon Foc inoculation. Time series profiling of proteome and 
metabolome of Foc-inoculated resistant and susceptible chickpea 
roots was performed using label-free quantitative proteomics and 
untargeted 'H-NMR metabolomics. We observed highly orches- 
trated response with significant modulation in various metabolic 
processes. The results described here thus improve our funda- 
mental knowledge of molecular dynamics associated with the 
chickpea-Foc interaction and potentially useful in designing 
strategies against wilt disease in chickpea. 


Results 


Protein identification and quantification in Foc- 
inoculated chickpea roots 


Two days after inoculation (DAI) with Foc, the susceptible 
chickpea cultivar (JG62) showed yellowing phenotype followed 
by drooping of leaves that finally lead to complete wilting by 12 
DAI. Whereas, Foc-inoculated resistant cultivar (Digvijay-DV) and 
the mock-inoculated wilt susceptible and resistant cultivars 
remained healthy throughout the experimental period. The 
JG62 plant could not sustain the fungal invasion beyond 12 
DAI, while DV plants remained unaffected (Figure S1). Based on 
these phenotypes, 2 and 4 DAI were considered as early stage, 
while 8 and 12 DAI were considered as late stage. We conducted 
high-throughput label-free quantitative proteomics analysis with 
Foc- and mock-inoculated chickpea root tissues at various time 
points from 2 to 12 DAI. This analysis identified a total of 811 
proteins (Table S1a) from which fungal proteins were excluded in 
further analysis (Lee et al., 2009). From these, 481 had statisti- 
cally significant differential expression (P < 0.05 and fold change 
>1.2) across cultivars and over the course of infection (Tables S1b 
and c). The ratio of normalized intensity of proteins from the 
inoculated samples vis-a-vis respective controls revealed increased 
or decreased expression in the Foc-inoculated roots. The log2- 
transformed values of differentially expressed proteins were 
clustered using SplineCluster (Heard et al., 2006), which gener- 
ated eight clusters using a prior precision of 1 x 1074 (Figure 1). 


Protein expression patterns in Foc-inoculated resistant 
and susceptible cultivars 


The potential biological function for each cluster was deduced 
based on gene ontology enrichment analysis using BiNGO. 
Cluster 1 (C1) had 61 proteins enriched for isoflavonoids 
biosynthesis and response to oxidative stress (Figure 1), while 
there were 73 proteins enriched in lignin biosynthesis, s-adenosyl 
methionine biosynthetic process and glycolytic process in C2. 
Majority of proteins from C1 and C2 revealed higher expression 
at all the stages except 8 DAI in DV compared to JG62. The C3 
cluster with 98 proteins was enriched for stress response, malate 
metabolism, oxidation—reduction process and glycolytic pathway. 
These proteins had increased expression at early stages which 


decreased at later stages in DV. Total 63 proteins from C4 cluster 
were enriched for the response to misfolded proteins, micro- 
tubule polymerization processes, osmotic stress response, pro- 
teosome core complex assembly and gluconeogenesis. The C5 
cluster had 44 proteins enriched for protein degradation through 
ubiquitin, ATP biosynthesis, photorespiration and active proton 
(H*) transport. Proteins from C4 and C5 clusters showed general 
trend of higher expression in all the stages except 12 DAI in DV. 
On the contrary, JG62 showed lower expression in majority of 
these proteins at later stages (Figure 1). Total 54 proteins had 
high expression in C6 at all the stages except 12 DAI in DV 
compared to JG62 and did not show enrichment to any specific 
process through BINGO. The C7 cluster had 52 proteins that were 
enriched in defence response and oxidative stress response. These 
proteins displayed overall contrasting expression patterns in DV 
and JG62 at early (2 DAI) and late (12 DAI) stages. The 36 proteins 
from C8 cluster were enriched in response to abiotic stress, fatty 
acid biosynthesis and nucleosome assembly processes. 


Foc induced quantitative variation in proteins from 
important metabolic pathways 


Upon Foc inoculation, we observed a complex response in 
chickpea from interconnected metabolic pathways including 
primary amino acid metabolism, glycolysis/gluconeogenesis, 
TCA cycle, phenylpropanoid pathway and increased lignification 
(Figure 2). Other cellular processes altered during Foc infection 
included unfolded protein response (UPR) and defence-related 
proteins. Enzymes such as, sucrose synthase, phosphoglucomu- 
tase, transaldolase, enolase, pyruvate dehydrogenase, citrate 
synthase, succinyl-coA ligase, fructose bis-phosphate aldolase, 
phosphogluco kinase, phosphoglyceratemutase, fumarate dehy- 
dratase and malate dehydrogenase were up-regulated up to 3.0- 
fold during early and late stages following Foc inoculation in DV. 
However, they showed up to 2.0-fold decrease in the roots of 
JG62. Other proteins from the same pathway such as fructok- 
inase, succinate dehydrogenase and aconitase showed up to 
2.56-fold higher level in DV, while JG62 exhibited reduction by 
1.4-fold. Several proteins from the UPR pathway including Hsp70, 
luminal binding protein (BiP), calmodulin, component of SCF-for 
SKP1-Cullin-F box protein and protein disulfide isomerase (PDI) 
increased up to 2.6-fold in DV, while they showed 2-fold 
reduction in JG62 upon Foc inoculation. The majority of these 
proteins were high in the early stage and low in the late stage in 
DV, while JG62 showed the opposite expression trend. Further, 
proteins induced by abiotic stress such as profilin and aquaporin 
PIP-type 7a revealed 1.7 to 2.5-fold increase in DV upon Foc 
inoculation compared with JG62. However, nuclear transcription 
factor-Y (NF-Y) had >30-fold increase at both the stages in JG62 
compared with DV. 

Proteomic analyses further depicted that proteins involved in 
reactive oxygen species (ROS) generation such as ascorbate 
peroxidase, peroxiredoxin, dehydroascorbate reductase (DHAR), 
hydroxyacyl glutathione hydrolase, glutathione peroxidase, 
glutaredoxin, glutathione S-transferase (GST), quinone oxidore- 
ductase and copper amine oxidase (CUAO) showed 3.0 to 6.0- 
fold increased expression in DV compared to JG62 exhibiting 
significant oxidative stress in chickpea after Foc inoculation. 
Additionally, enzymes involved in methionine metabolism, includ- 
ing methionine synthase, adenosine homocysteine hydrolase, 
adenosine kinase and AdoMet synthetase showed up to 1.6-fold 
higher levels in the Foc-inoculated DV at early and late stages, 
whereas JG62 displayed down-regulation by >1.3-fold at both 
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Figure 1 Clusters (C1 to C8) of 481 differential C6 (54) - 


expressed proteins in chickpea root. For each 
protein, the ratio of Logz-normalized expression 
of Foc inoculated with its respective control at 
various stages (2 to 12 DAI) and represented by a 
colour, according to the colour scale at the top. 
The number of proteins in a given cluster with 
similar expression trend is indicated in 
parentheses. v 
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the stages. In the present study, subsets of proteins participating 
in multiple branches of phenylpropanoid pathway such as lignin, 
flavonoid, isoflavonoid and phenolic biosynthesis were identified 
in Foc-inoculated chickpea cultivars. Caffeic acid O-methyltrans- 
ferase (CAOMT) and caffeoyl-CoA O-methyltransferase 
(CCoAMT) showed up to 2.0 and 3.85-fold increase, respectively, 
in DV compared to JG62. Similarly, enzymes from isoflavonoid 
biosynthesis such as chalcone synthase (CHS), chalcone isomerase 
(CHI), isoflavone synthase (IFS) and isoflavone reductase (IFR) 
revealed 1.8 to 5.0-fold increase in DV than JG62. 

Additionally, pathogenesis-related (PR) proteins such as endo 
B-1,3-glucanase, major latex protein (MLP), major latex allergen 
hev b5 and Bet v1 showed up to 2.1-fold higher expression in DV 
compared with JG62. Similarly, B-gulcosidase, disease resistance 
response (DRR) protein-206 and DRR-49 had >5.0-fold higher 
accumulation in DV than in JG62. Proteolytic chitinases offer 
antifungal properties and confer resistance to fungal pathogens. 
Chitinase, selenium binding protein (SBP) and glycine-rich pro- 
teins revealed up to 1.8-fold increased expression in DV 
compared with JG62. Interestingly, some PR proteins such as PR 
protein STH-2 (2-fold), thaumatin-like protein PR-5b (>5-fold) and 
PR-4A (>10-fold) were very high at early stage in DV than JG62. 
However, these proteins exhibited reverse trend with 2.5 to 10- 
fold increase in JG62 at late stages of Foc infection indicating 
their response to heavy wounding in chickpea. Some of the 
proteins involved in the stress signalling process such as ABA- 
responsive protein, auxin-binding protein ABP19a and Ran- 
binding protein were increased by 2.1-fold in DV compared with 
JG62. Additionally, 14-3-3 and H*-ATPase showed up to 2.0-fold 
higher levels in DV than in JG62, suggesting their role in defence 
response against Foc in chickpea roots. 
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Metabolic profiling in chickpea root 


A typical annotated 'H-NMR spectrum of chickpea root extract 
is depicted in Figure 3. The metabolite resonances were 
assigned according to the in-house databases and previous 
publication (Fan, 1996). These were further confirmed with a 
series of 2D NMR spectra including 'H-'H correlation spec- 
troscopy (COSY), 'H-'H total correlation spectroscopy (TOCSY), 
J-resolved spectroscopy (RES), 1H-13C heteronuclear single 
quantum coherence spectroscopy (HSQC) and 'H-'3C heteronu- 
clear multiple-bond correlation (HMBC) with both 'H and 1°C 
chemical shifts and signal multiplicities were as shown in 
Table S2. A total of 52 dominant metabolites were identified 
including amino acids (Ala, Val, lle, Leu, Gln, Glu, Asn, Trp, Lys 
and GABA), sugars (glucose, sucrose, fructose, trehalose and 
salicin), organic acids (pyruvate, lactate, acetate, citrate, succi- 
nate, formate, fumarate, malate and guanidoacetate), nucleo- 
sides (adenosine, uridine, inosine, 5CMP and hypoxanthine) and 
phytoalexins (genistein and luteolin). The identification of few 
of these metabolites was confirmed by spiking with their 
known standards (Table S2). 


Differential metabolic alterations induced by Foc 
inoculation 


Principal component analysis (PCA) revealed significant metabolic 
changes between Foc-inoculated resistant and susceptible culti- 
vars at different time points. The averaged PCA scores were 
calculated for the first two PCs to construct the PCA trajectory that 
illustrated clear separation in metabolites from control and 
inoculated samples between two cultivars (Figure 4). Interestingly, 
the metabolic profiles obtained from the controls of both cultivars 
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Figure 2 Interconnection between various metabolic processes in chickpea—Foc interaction. Each graph represents differential expression (fold 
change) pattern according to the colour scale; columns represent the four stages (2 to 12 DAI) after Foc inoculation. The first row represents the 


resistant cultivar, while the second row represents the susceptible cultivar. 


and the inoculated DV followed similar trajectory trends; however, 
the inoculated JG62 showed dramatic change in trajectory after 8 
DAI (Figure 4). The metabolic changes resulted after Foc inocu- 
lations were further evaluated by constructing orthogonal projec- 
tion to latent structures discriminant analysis (OPLS-DA) models. 
The quality of models was indicated by the values of R? and Q? and 
cross-validated with a CV-ANOVA approach (P < 0.05) (Figures 5a 
and b) and permutation tests (Figure S2). All the significantly 
changed metabolites are annotated in the OPLS-DA coefficient 
plots (Figures 5a and b) and summarized in Table $3. 

Compared with the respective control, Foc inoculation induced 
marked reduction in the levels of amino acids including Thr, Ala, 
Lys and Asn at early stage in both the cultivars. However, their 
levels were significantly elevated in JG62 at late stage while DV 
remained unchanged (Figure 5). The decreased levels of Ile and 
Leu were observed at early stage and remained unchanged at 
late stage in JG62 after Foc inoculation. Interestingly, Val and Trp 
levels increased dramatically in the late stage of JG62 but 
remained unchanged in DV at both the stages (Figure 5). 
Similarly, two of the most important metabolites from nitrogen 
metabolism, such as Glu and Gln, were up-regulated only at early 


stage of Foc inoculation in DV. However, JG62 revealed 
significantly higher accumulation of both of these amino acids 
at early and late stages. In addition, Foc inoculation resulted in 
decreased level of glucose in DV at early stage; however, no 
significant change was observed in both the cultivars at late 
stage. The levels of sucrose and fructose were decreased at early 
stage in JG62, followed by significant reduction at the later stage. 
t is of particular interest that the levels of some nucleotides 
including uridine and orotate were significantly decreased at the 
ater stage of JG62 compared with DV. However, the opposite 
trend in the levels of adenosine and inosine was observed at the 
ater stage of JG62. Compared with the controls, Foc inoculation 
induced significant reduction in the levels of malate and acetate 
in JG62 at the later stage. As for phytoalexins, genistein was 
increased in DV at the later stage but remained unchanged in 
JG62 and the level of luteolin decreased at the later stage of 
JG62 while DV was unaffected. A known antifungal compound, 
clotrimazole level increased in DV, but there was a marked 
reduction in JG62 at the later stage. In addition, Foc infection 
caused increase in quinone at early stage and phytosterol at the 
later stage of DV only. 
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Figure 3 Typical NMR spectrum of chickpea root extract. Annotation with number and details of metabolites is provided in Table S2. 
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Figure 4 PCA trajectory plots. (a) resistant-DV and (b) susceptible-JG62 plants with their respective controls obtained from mean of PC1 and PC2 
values at 2 to 12 DAI with error bars representing two standard deviations. Foc-inoculated samples are in green while respective controls in red. Top 


right corner box indicates overall pattern. 


Comparative expression of candidate genes in root 


To obtain complementary information of transcriptional varia- 
tions, we examined expression levels of key genes (Figure 6) 
involved in various metabolic pathways based on our proteomic 
and metabolomic findings. These included genes from nitrogen 


mobilization (glutamate dehydrogenase-GDH, glutamate synthase, 
glutamine synthase and asparagine synthase), stress response 
(NF-Y and SKP1-like protein 1A), methionine metabolism 
(methionine synthase and AdoMet synthetase), lignin and 
phytoalexin biosynthetic pathways (CCoAMT, CHS, CHI, iso- 
flavone 4’-O-methyltransferase, IFS and IFR). GDH expression 
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Figure 5 Pairwise comparison via OPLS-DA. OPLS-DA scores plots (left) and corresponding coefficient-coded loadings plots (right) obtained from 
metabolic profiles of Foc-inoculated (a-d) resistant-DV and (e-h) susceptible-JG62 cultivars and their respective controls at 2 to 12 DAI. The coloured scale 
in correlation coefficient (|r|: absolute values) plots shows the significance of metabolite variations discriminating between the Foc-inoculated and control 


plants. 


markedly increased by >200 fold in the inoculated JG62 at the 


later stage. Similarly, significant up-regulation (up to 4-fold) of 
glutamine synthetase, asparagine synthetase and glutamate 
synthase was observed in the Foc-inoculated JG62 (Figure 6a— 
d). However, candidate genes from methionine metabolism 
such as Adomet synthetase and methionine synthase revealed 
>3 and 10-fold enhanced expression, respectively, in DV 
(Figure 6e and f). Also, important genes from lignin and 
phytoalexin biosynthetic pathway displayed higher expression in 
DV compared with JG62 (Figures 6i-n). On the contrary, abiotic 
stress induced NF-Y gene showed >8-fold higher expression in 


JG62 (Figure 6h), 
(Figure 6g). 


while SKP1 was expressed more in DV 


Intense lignification is associated with Foc resistant 
phenotype 


After Foc inoculation, differential lignin deposition in the root 
tissue of DV and JG62 was observed as the time progressed, 
wherein DV exhibited intense lignification compared with JG62 
(Figure S3). The fungal pathogen invaded the susceptible cultivar 
and blocked the vascular tissue completely by 12 DAI, leading to 
wilting of JG62 plants. 
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Figure 5 Continued. 


Discussion 


Foc induced remodelling in energy metabolism and 


nitrogen mobilization 


Obligate biotrophs depend on host metabolism for nutrient uptake, 
which in turn determine their pathogenicity within the host. In the 
present study, proteins involved in glycolysis and TCA cycle were 


Foc infection in both the cultivars; however, 


this was more 


predominant in the JG62. Similar rapid reduction in the levels of 


these sugars was also reported in sunflower upon infection with 


Botrytis cinerea (Dulermo et al., 2009b). These results emphasized 


up-regulated in DV while down-regulated in JG62. We also 
observed steep alteration in primary metabolites (amino acids and 


sugars) specifically in JG62. Sugars affect disease susceptibility 
often favouring disease development while playing a critical role in 


innate defence pathways involving metabolic regulation (Bolouri 
Moghaddam and Van den Ende, 2012). In the present study, rapid 


decrease in sugars such as sucrose and fructose was observed upon 


the regulatory role of sugars in the metabolic reprogramming and 
subsequently higher expression of proteins from glycolysis and TCA 
cycle endowed DV plants to successfully combat Foc invasion. 
Moreover, supportive evidence could be derived from earlier 
reports of F. oxysporum induced up-regulation of various ESTs 
from sugar metabolism in chickpea (Ashraf et al., 2009; Gupta 
et al., 2010, 2013a). 


Nitrogen plays an essential role in the nutrient relationship 


between plants and pathogens. 
concentration due to nitrogen mobilization after pathogen 


Modulation in amino acid 
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Figure 6 Quantitative real-time PCR of various candidate genes. (a-n) expression variation in each gene observed from chickpea root of 
Foc-inoculated resistant-DV and susceptible-JG62 cultivars as compared to their respective controls at 2 to 12 DAI. 


© 2016 The Authors. Plant Biotechnology Journal published by Society for Experimental Biology and The Association of Applied Biologists and John Wiley & Sons Ltd., 1- 


(i) Chalcone synthase 


(k) Chalcone isomerase 


Relative transcript abundance 
on 
Gr 
Coy 
“Cry 


1 I I 
i<-- Susceptible --> |<- Resistant ->! 
! ! ! 


(m) Caffeoyl-CoA O-methyltransferase 
0.05 


0.04 
0.03 
0.02 


0.01 Al 
= a —= 


0,00035 
0.00030 
0.00025 
0.00020 
0.00015 
0.00010 
0.00005 
0.00000 


Relative transcript abundance 


Susceptible --» aks 


Figure 6 Continued. 


infection to plant has been reported (Dulermo et al., 2009a; 
Tavernier et al., 2007). Proteomics studies on wheat-F. gramin- 
earum exhibited increase expression in proteins from amino acid, 
carbon and nitrogen metabolism (Wang et al., 2005; Zhou et al., 
2006). Such alteration with significantly low amino acid content 
was also observed in tomato and sunflower after infection of 
B. cinerea and Sclerotinia sclerotiorum, respectively (Berger et al., 
2004; Jobic et al., 2007). Current study showed significant 
decrease in the concentration of various amino acids in the 
susceptible cultivar, which suggested that the fungus probably 
utilized them for its establishment and proliferation inside the 
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host. However, probably due to Foc sporulation, levels of amino 
acids were increased at the later stage in JG62. Hence, the role of 
nitrogen mobilization was further scrutinized in chickpea—Foc 
interaction by gene-expression analysis of four representative 
enzymes from nitrogen mobilization viz. glutamate synthase, 
glutamate dehydrogenase, glutamine synthetase and asparagine 
synthetase. Significant up-regulation of these genes in JG62 
compared to DV correlated well with the metabolomics outcome. 
Consistently, proteomic analysis also revealed higher level of 
sucrose synthase and glutamine synthase following Foc inocula- 
tion. Thus, proteomic, metabolomic and gene-expression results 
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together indicated that the process of nitrogen mobilization could 
be critical in the establishment of Foc infection in the susceptible 
plants. 


Stress responsive proteins in chickpea root 


During stress, accumulation of unfolded proteins increases in 
endoplasmic reticulum (ER) and results in triggering the 
unfolded protein response (UPR) to remove the misfolded 
proteins by ubiquitin-proteasome pathway. Thus, UPR not only 
helps to avert the cytotoxic impact of misfolded proteins, but 
also assists to relieve stress and reinstate normal functions in ER 
(Ye etal., 2011). Up-regulation of critical proteins from UPR 
pathway such as Hsp70, BiP, calmodulin, SKP1 and PDI in DV 
could suggest their coordinate response. In Arabidopsis 
thaliana roots, swelling of ER and vacuolar collapse resulting 
in ER stress and cell death were observed during fungal 
colonization (Qiang et a/., 2012). Similarly, PDI level was higher 
in wheat plants upon Puccinia striiformis inoculation (Maytalman 
et al., 2013). This study also reinforced that during defence 
response, there was constant requirement for proteins stabiliza- 
tion in the process of folding, assembly, vesicle trafficking and 
secretion. Collectively, these outcomes suggested the impor- 
tance of an efficient UPR pathway utilization in the resistant 
chickpea against Foc. 

During stress conditions, aquaporins such as PIP are involved 
in water transport in plant. Current investigation revealed 
significant increase in PIP-7a level at both the stages in DV, 
while drastic reduction was observed in JG62. This could have 
resulted in better water conductance in DV and helped the 
resistant plant to overcome fungal attack. On the contrary, 
diminished water transport with reduced aquaporin in xylem 
due to fungal spread in the susceptible plant could have 
eventually turned into wilting. Furthermore under water-limited 
conditions, nuclear factor-Y (NF-Y) family and ABA-responsive 
protein have shown to be up-regulated in Arabidopsis (Nelson 
et al., 2007). We observed high expression of these proteins at 
both the stages in JG62 compared with DV. Thus, present 
findings indicated constrain in water uptake in JG62 due to 
clogged xylem after fungal invasion that potentially increased 
the susceptibility to Foc. 


Early recognition of Foc leads to ROS generation and 
lignosuberization 


Generation of ROS is one of the earliest cellular responses to 
pathogen recognition and/or infection (Gupta et al., 2013b). The 
enzymes involved in ROS production such as peroxidase, DHAR, 
hydroxyacyl glutathione hydrolase, glutathione peroxidase, 
glutaredoxin, GST, quinone oxidoreductase and CuAO were 
significantly increased in DV compared with JG62 in the current 
investigation. It has been demonstrated that CuAO and perox- 
idases functionally correlate in lignosuberization process (Angelini 
et al., 1993; Scalet et al., 1991) and the inhibitors of CuAO result 
in decreased defence response (Rea et al., 1998, 2002). Addi- 
tionally, Raju et al. (2008) found more lignification in the resistant 
chickpea cultivar compared with the susceptible one upon Foc 
infection. In the present study also, intense lignin deposition on 
the root cortex of Foc-inoculated resistant cultivar was observed 
(Figure $3). Taken together, ROS generation and higher expres- 
sion of CuAO suggested that Foc triggered hydrogen peroxide 
generation and lignosuberization process leading to initiation of 
defence response in DV. Secondly, monolignol biosynthesis also 
plays critical role in the host defence mechanism through 


lignification, making cell wall more resistant to the pathogen 
penetration. Reduction in the monolignol biosynthesis through 
co-silencing of CAOMT and CCoAMT in wheat led to the higher 
penetration by Blumeria graminis (Bhuiyan et al., 2009). We also 
found up-regulation of these two lignin biosynthetic enzymes in 
DV than JG62 suggesting increased lignin deposition that might 
provide resistance against Foc. 


Crucial role of methionine metabolism in Foc resistance 


In the Yang cycle, methionine synthase converts homocysteine to 
methionine contributing in synthesis of Adomet by AdoMet 
synthetase. AdoMet can also lead to ethylene by 1-aminocyclo- 
propane-1-carboxylic acid (ACC) synthase and ACC oxidase. 
However, down-regulation of ACC oxidase by 1.5-fold in both 
the cultivars indicated that preferential AdoMet pool was not 
channelized towards ethylene production. A recent proteomic 
study has shown that increased expression of AdoMet synthetase 
and lower level of ACC oxidase resulted in higher methionine 
recycling and lower ethylene biosynthesis in rice roots against 
Herbaspirillum seropedicae (Alberton et al., 2013). Silencing of 
AdoHcy hydrolase in transgenic tobacco plants confirmed its role 
in defence mechanism against pathogens (Masuta et al., 1995). 
Similarly, Kawalleck et al. (1992) identified mRNAs for AdoMet 
synthetase and AdoHcy hydrolase in parsley plant upon fungal 
infection. Altogether, current study indicated close association 
between pathogen defence and increased level of activated 
methyl groups during chickpea—Foc interplay. 

Further, highly methylesterified pectin is required for normal 
plant cell wall, while it is de-esterified by pectin methyl esterase 
(PME) which leads to increased vulnerability of plant cell wall to 
the pathogen invasion (Lionetti et al., 2007). Our proteomics 
results showed increased expression of PME at both the stages in 
JG62 compared with DV. Earlier studies with either silencing of 
PME or overexpression of PME inhibitors in plants demonstrated 
negative role of PME in the pathogen resistance (An et al., 2008; 
Ma et al., 2013). Likewise, plant sterols are structurally related to 
cholesterol and control mechanical property of cell membrane 
(Hodzic et al., 2008) and also serve as substrate for several 
metabolic pathways (Itkin etal, 2013). Interestingly, our 
metabolomic data also revealed high phytosterol at later stages 
in DV. However, there was no significant change in JG62 except 
at 8 DAI. The stability of the plant cell wall is also maintained by 
actin binding cytoskeleton protein such as Profilin, which showed 
higher expression in DV compared with JG62 at late stage. 
Overall, altered methionine metabolism affected methyl esterifi- 
cation of pectin in Foc-infected susceptible roots, while stronger 
plant cell wall with normal pectin and cytoskeleton proteins could 
have helped the resistant chickpea cultivar with better defence 
against Foc. 


Foc resistance is mediated by phenylpropanoid pathway 


Plants respond to pathogen challenge by increased activation of 
the phenylpropanoid pathway leading to flavanoids, isoflavo- 
noids and phenolics biosynthesis. They play multiple roles in 
plant-pathogen interaction including precursors for the defence- 
related phytoalexins and signal molecules in response to 
pathogen infection. As detailed in the results, enzymes involved 
in this pathway such as CHS, CHI, IFR and IFS were up-regulated 
in DV as compared to JG62. In the transgenic soybean roots, 
RNAi silencing of CHS gene showed decrease in total isoflavo- 
noids as well as reduced resistance to fungal pathogens 
(Lozovaya et al., 2006; Subramanian et al., 2005). Additionally, 
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Naoumkina et al. (2007) and Farag et al. (2008) showed that 
fungal extract induced higher levels of CHI, IFS and IFR leading to 
the production of phytoalexin in Medicago cell suspension 
culture to combat the pathogen. Consistently, the accumulation 
of genistein, luteolin and quinone identified by metabolome 
analysis in the present study correlated well with the proteomic 
data. Thus, the accumulation of isoflavonoid biosynthetic pro- 
teins and metabolites in DV suggested their potential involve- 
ment in Foc resistance. 


Modulation of defence-related proteins and metabolites 
upon Foc inoculation 


In the present study, we observed quantitative variation in many 
defence-related proteins such as endo §-1,3-glucanase, MLP, 
hev b5 and Bet v1, B-gulcosidase, DRR-206, DRR-49, chtinases, 
SBP, PR10, STH-2, PR4a, PR-5b, 14-3-3 and H*-ATPase in Foc- 
inoculated chickpea cultivars. Many previous studies in other 
plant-pathogen interactions have also demonstrated the impor- 
tance of such proteins in plant defence. Lytle et al. (2009) and 
Gurjar etal. (2012) have reported that Bet v1 protein is 
responsible for resistance towards pathogen. Similarly, B-1,3- 
glucanase is involved in plant defence response against 
pathogen (Shetty et al., 2009; Ward et al., 1991). Previously 
two independent studies have shown that B-glucosidase, DRR- 
206 and DRR-49 proteins contribute to lignification process 
(Burlat et al., 2001; Hosel and Barz, 1975). Increased chitinase 
activity was also reported earlier after Foc inoculation in the 
resistant chickpea cultivar (Giri et al., 1998). Similarly, overex- 
pressed SBP in rice provided more resistance against rice blast 
fungus (Sawada et al., 2004). Another important defence- 
related protein, 14-3-3, is known to be associated with 
hypersensitive cell death in pathogen incompatible cultivars 
(Roberts, 2003). In our earlier study, up-regulation of 14-3-3 
transcripts was detected in Foc-inoculated resistant chickpea 
roots (Nimbalkar et al., 2006). Thus, up-regulation of 14-3-3 
and H*-ATPase suggested their roles in activating hypersensi- 
tivity response in chickpea leading to Foc resistance. Apart from 
the above-mentioned proteins, our metabolite analysis revealed 
increase level of clotrimazole in Foc-inoculated DV at late stage 
compared to JG62. This is an interesting observation, clotrima- 
zole being an anthropogenic antifungal compound and needs 
further experimentation. Thus, the pathogen attack triggered 
expression of defence-related genes, secondary metabolites 
with antimicrobial nature and PR proteins in the chickpea—Foc 
interactions. 

In summary, our integrated approach of label-free quantita- 
tive proteomics, 'H-NMR metabolomics and candidate gene- 
expression analysis has facilitated to understand the defence 
mechanism in chickpea against Foc. Most of the proteins/ 
metabolites from various metabolic processes such as energy 
metabolism, isoflavonoid/flavonoid biosynthesis pathway and 
lignin biosynthesis that lead to defence were up-regulated 
in the resistant plant compared with the susceptible one. 
Further, elevated levels of some of the proteins from the UPR 
pathway during incompatible interaction indicated that proper 
protein folding and transport might be crucial for the plant's 
survival. Some of the proteins/metabolites such as isoflavone 
reductase and isoflavone synthase leading increased phytoalex- 
ins levels, CAOMT and CCoMT involved in lignosuberization 
and increased methionine synthase for efficient methylation 
process in the resistant plant could have helped against the 
pathogen. Overall, above findings improve our understanding 
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on the metabolic reprogramming during the wilt disease 
progression. 


Experimental procedures 
Plant material 


Seeds of wilt resistant (Digvijay, DV) and susceptible JG62, JG) 
chickpea cultivars as well as pathogenic culture of Foc race 1 were 
obtained from Mahatma Phule Krishi Vidyapeeth (MPKV), Rahuri, 
India. All the steps of plant growth and pathogen inoculation 
were followed as previously described (Kumar et al., 2015) and 
divided into two groups, control or mock (sterile water) inocu- 
lated and Foc (~10® spores/mL) inoculated. The workflow for 
sample collection/data analysis is shown in Figure S4. As the 
pathogen is reported to colonize the xylem vessels 2 days 
postinoculation (Gupta etal., 2010), the root tissues were 
collected at 2, 4, 8 and 12 DAI. Whole roots of chickpea was 
used in equal amount from the resistant and susceptible cultivars 
at specific time points simultaneously for the proteomics and 
metabolomics analysis; because the fungal hyphae might not be 
localized but distributed in the whole root xylem tissue (Jimenez- 
Fernandez et a/., 2013). In each group, a pool of 10 plants 
comprised one biological replicate. For proteomics analyses, 
tissues from three biological replicates were analysed, while ten 
biological replicates were analysed for metabolomics analysis. The 
harvested tissues were stored at —80 °C till further use. 


Protein extraction and mass spectrometry analysis 


Total proteome of root tissue was extracted as described by 
Isaacson et al. (2006). Protein pellets were solubilized in 50 mm 
ammonium bicarbonate buffer containing 0.1% Rapigest 
(Waters, Milford, MA). The dissolved proteins were reduced and 
alkylated by DTT and iodoacetamide, respectively, followed by 
overnight tryptic hydrolysis at 37 °C using Promega sequencing 
grade trypsin. The digested peptides were analysed with LC-MSE 
workflow using nano-ACQUITY online coupled to a SYNAPT 
HDMS system (Waters). Nano-LC separation was performed with 
symmetry C18 trapping column (180 um x 20 mm, 5 um) and 
bridged-ethyl hybrid (BEH) C18 analytical column 
(75 um x 250 mm, 1.7 pm). The binary solvent system com- 
prised solvent A (0.1% formic acid in water), and solvent B (0.1% 
formic acid in acetonitrile). Each sample (500 ng) was initially 
applied to the trapping column and desalted by flushing with 1% 
solvent B for 1 min at a flow rate of 15 L/min. Elution of the 
tryptic digested sample was performed at a flow rate of 300 nL/ 
min by increasing the solvent B concentration from 3% to 40% 
over 90 min. Before data acquisition, the mass analyser was 
calibrated using Glu-fibrinopeptide B (Sigma-Aldrich, Steinheim, 
Germany) from m/z 50 to 1990. The Glu-fibrinopeptide B (GFP-B) 
was delivered at 500 fmole/uL to the mass spectrometer via a 
NanoLockSpray interface using the auxiliary pump of the nano- 
ACQUITY system at every 30 s interval for lock mass correction 
during data acquisition. Data-independent acquisition was per- 
formed (LC-MS) as described by Patel et al. (2009). 

As accuracy and reproducibility in mass measurement are 
critical in data acquisition during large-scale proteomic experi- 
ments, PCA was used to assess the quality of the measurement 
in terms of replicate similarity. The replicates of each sample 
were clustered together reflecting inherent similarities between 
samples (Figure S5a). In addition, linear response and repro- 
ducibility of measurement of the quantitative proteomic data 
acquisition were tested by plotting two replicates (Figure S5b), 
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whereas Figure S5c showed that data have been acquired below 
3 ppm mass accuracy. Further, the percent coefficient of 
variance of retention time (% CV-RT) was calculated to assess 
the separation stability and coefficient of variance of 0.3 min, 
which also suggested stability in chromatographic separation 
(Figure S5d). 


Analysis of quantitative proteomics data 


The acquired LC-MS® data were processed using the ProgenesisQI 
for Proteomics software (Waters). Protein identifications were 
obtained by searching the genomic databases of chickpea (http:// 
www. icrisat.org/) and Fusarium oxysporum (http://www.broad- 
institute.org/). LC-MS® data were searched with a fixed car- 
bamidomethyl modification for cysteine residues, along with a 
variable modification for oxidation of methionine, N-terminal 
acetylation, deamination of asparagine and glutamine and 
phosphorylation of serine, threonine and tyrosine. The ion 
accounting search algorithm within ProgenesisQI for Proteomics 
software was used which has been developed specifically for 
searching data-independent MS® data sets and described by Li 
et al. (2009). The ion accounting search parameters were (a) 
precursor and product ion tolerance: automatic setting, (b) 
minimum number of product ion matches per peptide: 3, (c) 
minimum number of product ion matches per protein: 7, (d) 
minimum number of peptide matches per protein: 1 and (e) 
missed tryptic cleavage sites: 1. False-positive rate was set at 1%. 
Search results of the proteins and the individual MS/MS spectra 
with a confidence level at or >95% were accepted. Label-free 
quantitation of identified protein was done on the basis of spiked 
bovine serum albumin (BSA) protein. 


Clustering of identified proteins 


Data were normalized by spiked BSA (50 fmoles), and relative 
accumulation differences were determined for proteins having 
differential expression. Sum of three replicates of inoculated 
samples was divided by that of the respective controls. This 
established a ratio of accumulation of a protein in plants upon Foc 
inoculation compared with mock-inoculated control plants. The 
log2-transformed ratio (susceptible/control and resistant/control) 
pairs were clustered by the application of SplineCluster (Heard 
et al., 2006), a Bayesian model-based hierarchical clustering 
algorithm for time series data. 


Gene ontology enrichment analysis 


Protein functional annotation was determined using Blast2GO 
(Conesa etal., 2005) and for each cluster, GO enrichment 
analysis was carried out using BINGO 2.3 plugin tool in Cytoscape 
version 2.8 (Maere et al., 2005). Over-represented GO biological 
process categories were identified using a hypergeometric test 
with a significance threshold of 0.05 after Benjamini and 
Hochberg false discovery rate correction (Benjamini and Hoch- 
berg, 1995) using the annotated chickpea genome as the 
reference set. 


Metabolite extraction and NMR measurement 


Plant root tissue was ground well in liquid nitrogen by using 
bead beater (Retsch GmbH, Retsch-Allee, Germany) and 
lyophilized. The powdered root tissues (~50 mg) were extracted 
with 0.75 mL of CD30D and 0.75 mL of 10 mm KH>PO, buffer 
(pH 6.0) containing sodium3-trimethlysilyl [2,2,3,3-D,4] propi- 
onate (TSP) as described previously (Kim et al., 2010). After 
ultrasonication for 20 min and centrifugation at 12 000 g for 


10 min at room temperature (~25 °C), 0.5 mL of supernatant 
was collected for NMR detection. 'H NMR spectra of root 
extracts were acquired at 25°C on a Bruker AV Il 500 
spectrometer (Bruker Biospin, Rheinstetten, Germany) operating 
at 500.13 MHz for 'H. A standard water-suppressed one- 
dimensional NMR spectrum was recorded using noesyprid 
pulse sequence (RD-90°-t;-90°-tm-90°-acquisition) with the 
recycle delay of 6s and the mixing time (tm), of 50 ms. 
Typically, 90° pulse was set to about 15 us and 256 transients 
were collected into 48K data points for each spectrum with a 
spectral width of 16 ppm. All spectra were referenced to 
chemical shift of TSP (6 = 0.00). For the metabolite assignment 
purpose, a range of two-dimensional NMR spectra were 
recorded for selected samples including COSY, TOCSY, HSQC 
and HMBC. In COSY and TOCSY experiments, respective 64 
and 32 transients were collected into 2K data points for each 
of 256 increments with the spectral width of 2426 Hz for both 
dimensions. Magnitude mode was used with gradient selection 
for the COSY experiments, whereas the mlevgpphw5 pulse 
program was employed as the spin-lock scheme in the phase 
sensitive mode, with the mixing time of 60 ms, for TOCSY. 
Both HSQC and HMBC spectra were acquired using the 
gradient-selected sequences. In HSQC experiment, 80 transients 
were collected into 1K data points for each of 140 increments. 
In HMBC experiment, 160 transients were collected into 2K 
data points for each of 256 increments. The spectral widths 
were 2426 Hz for 'H and 9809 Hz for '?C in HSQC and HMBC 
experiments. Confirmation of resonance assignments were 
performed for few metabolites by spiking the samples with 
known standards (Table $2). 


NMR spectra processing and multivariate data analysis 


All the 'H NMR spectra were manually corrected for phase and 
baseline distortions using TOPSPIN (v2.1; Bruker Biospin) and 
calibrated for chemical shift drifting by in-house-developed 
script for MATLAB (The Mathworks, Natick, MA). The spectral 
region 6 0.5-9.5 was divided into bins with width of 0.002 ppm 
(1.0 Hz) using AMIX software (v3.8.3; Bruker Biospin GmbH, 
Germany). The region 6 4.727-5.089 ppm was discarded to 
remove the effects of imperfect water presaturation. The areas 
of the remaining bins were normalized to total sum of intensity 
for each spectrum to compensate for the overall concentration 
differences prior to statistical data analysis. Multivariate data 
analyses were carried out with SIMCA-P+ v 12.0 software 
package (Umetrics, Umea, Sweden). PCA was performed on the 
mean-centred NMR data to inspect overall data distributions and 
possible outliers. Using the NMR data as the X-matrix and group 
information as Y-matrix, OPLS-DA was carried out with unit 
variance scaling (Trygg, 2002; Xiao et al., 2008). The OPLS-DA 
models were 7-fold cross-validated and the quality of the model 
was described by the parameters R*X, representing the total 
explained metabolic variables and Q?, indicating the model 
predictability. The models were further evaluated with a CV- 
ANOVA approach (P < 0.05) and permutation tests. To facilitate 
interpretation of the results, back-transformation (Cloarec et al., 
2005) of the loadings generated from the OPLS-DA was 
performed prior to generating the loadings plots, which were 
colour-coded with the Pearson linear correlation coefficients of 
variables (or metabolites) using an in-house-developed script for 
MATLAB (The Mathworks) (Wang et al., 2007). The colour- 
coded correlation coefficient indicates the significance of the 
metabolite contribution to the class separation, with hot colours 
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(e.g. red) being more significance of the metabolite contribu- 
tions to the group classification than cold ones (e.g. blue). In this 
study, a correlation coefficient cut-off value of 0.602 (i.e., 
N= 10, |r| > 0.602) was used for the statistical significance 
based on the discrimination significance at the level of P < 0.05, 
which was determined according to the discriminating signifi- 
cance of the Pearson’s product-moment correlation coefficient 
(Cloarec et al., 2005). 


Quantitative real-time PCR analysis 


Total RNA was extracted from 100 mg root tissue by using TRI 
Reagent (Sigma-Aldrich). First strand cDNA synthesis was per- 
formed using the High Capacity cDNA Reverse Transcription Kits 
(Applied Biosystems, Foster City, CA) with 3 ug of DNase1 
treated total RNA using oligo (dT) primer following manufac- 
turer's protocol. Gene-specific primers were designed using 
Primer Express (v2.0) software (Applied Biosystems) and listed in 
Table S4. Real-time PCR was carried out as earlier described 
(Barvkar et al., 2012) using FastStart universal SYBR green master 
mix (Roche, Mannheim, Germany) with 7900HT Fast real-time 
PCR system (Applied Biosystems). The Initiation factor 4a (/F4cx) 
gene was used as internal standard or reference gene (Garg 
et al., 2010). 


Lignin staining 


Lignin accumulation within root tissue after Foc inoculation was 
detected using the phloroglucinol/hydrochloric acid stain as 
described by MauchMani and Slusarenko (1996). The control 
and Foc-inoculated chickpea roots were subjected to transverse 
sections with a scalpel and immersed in 1 mL of 1% phloroglu- 
cinol in 6N HCI for 5 min, and the lignin staining was visualized 
under light microscope. 
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